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ABSTRACT 

Nuclear rings in barred galaxies are sites of active star formation. We use hydrodynamic simulations 
to study temporal and spatial behavior of star formation occurring in nuclear rings of barred galaxies 
where radial gas inflows are triggered solely by a bar potential. The star formation recipes include a 
density threshold, an efficiency, conversion of gas to star particles, and delayed momentum feedback 
via supernova explosions. We find that star formation rate (SFR) in a nuclear ring is roughly equal to 
the mass inflow rate to the ring, while it has a weak dependence on the total gas mass in the ring. The 
SFR typically exhibits a strong primary burst followed by weak secondary bursts before declining to 
very small values. The primary burst is associated with the rapid gas infall to the ring due to the bar 
growth, while the secondary bursts are caused by re-infall of the ejected gas from the primary burst. 
While star formation in observed rings persists episodically over a few Gyr, the duration of active star 
formation in our models lasts for only about a half of the bar growth time, suggesting that the bar 
potential alone is unlikely responsible for gas supply to the rings. When the SFR is low, most star 
formation occurs at the contact points between the ring and the dust lanes, leading to an azimuthal 
age gradient of young star clusters. When the SFR is large, on the other hand, star formation is 
randomly distributed over the whole circumference of the ring, resulting in no apparent azimuthal age 
gradient. Since the ring shrinks in size with time, star clusters also exhibit a radial age gradient, with 
younger clusters found closer to the ring. The cluster mass function is well described by a power law, 
with a slope depending on the SFR. Giant gas clouds in the rings have supersonic internal velocity 
dispersions and are gravitationally bound. 

Subject headings: galaxies: ISM — galaxies: kinematics and dynamics — galaxies: nuclei — galaxies: 
spiral — ISM: general — shock waves — stars: formation 
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1. INTRODUCTION 



Nuclear rings in barred galaxies are sites of intense 
star formation (e.g. , |Burbidge fc Burbi dge|1960[| Sandage 
"^^^^^ |Phillips|1996| IButa & Combes"l996 Knapen et a' 



1961 



2006 



IMazzuca et al.||2008i iComeron et al.i|2010i |Sand- 



al.|2010[|Mazzuca et al.|2011[|Hsieh et al.|201ip : 
igs are tnought to form as a result ot nonlinear 



strom et al 
I'hese ringi 
interactio ns of gas with a non-axisymmetric bar p otcn- 



tial 



e^ 



1990 



Combes fc Gerin 1985; Shlos man et al.. 
Athana ssoulalll992j |Heller fc Shlosman 1994( IKn apen et 
al 111995; Buta & Combes 1996; Combes 2001; Finer et 



ar]|1995f [Regan &: Teuben|2003| ). Due to the bar torque, 
tHe gas readily torms dust-lane shocks in the bar region 
and flows inward along the dust lanes. The inflowing 
gas speeds up gradually in the azimuthal direction as it 
moves inward, and shapes into a ring very close to the 
galaxy center (e.g., Kim et al. |2012a ). Consequently, 
nuclear rings have very large surface densities and short 
dynamical time scales, capable of triggering starburst ac- 
tivity. 

There are some important observational results that 
may provide clues as to how star formation occurs in the 
nuclear rings. First of all, observations indicate that the 
star formation rate (SFR) in the nuclear rings appears to 
vary with time, and differs considerably from galaxy to 
galaxy. Analyses of various population synthesis mod- 

|seowy@astro.snu.ac.kr, wkim@astro.snu.ac.kr| 



els for a sample of galaxies reveal that the strength of 
observed emission lines from the nuclear rings is best 
described by multiple starburst activities over the last 
0.5 Gyr or so , rather than by a constant star fo rma- 
tion rate (e.g., |Allard et al.||2006| |Sg£gi et al.||2007[ ). For 
22 nuclear ringsT'Maz zuca et al.| ( 2D08| ) touiid that the 
present SFRs are widely distributed in the range 0.1- 
10 M0 yr~^. It appears that the SFR is largely insen- 
sitive to the total molecular mass in a ring, but can be 
strongly affected by the bar strength. For instance, the 
ring in a strongly-barred galaxy NGC 4314 has a rela- 
tively low SFR at ~ 0.1 Mq yr-i ( [Benedict et al.|2002J , 
which is about an order of magnitude smaller than tha t 



in a weakly-barred galaxy NGC 1326 ( |Buta et aIT2000[ ) 
although the total molecular mass contained in the ring 
is within a factor of two. In fact, the SFRs given in 
Mazzuca et aL (12008) combiiied with the bar strength 



presented by Comeron et al. (2010) show that strongly 
barred galaxies tend to have a very small SFR in the 
rings, while weakly-barred galaxies have a wide range of 
the SFRs. 

Second, based on the spatial distributions o f star- 
forming regions in nuclear rings, Boker et al. (20081 



"popcorn" and 



proposed two models of star formation: 

"pear ls on a string" models (see also Sandstrom et al. 

2010|. In the first popcorn model, star formation oc- 



curs m dense clumps that are randomly distributed along 
a nuclear ring. This type of star formation, presum- 
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ably caused by gr avitational instability of the ring itself 
( Elinegreen||1994 ), does not produce a systematic gradi- 
ent in the ages ot young star clusters along the azimuthal 
dir ection (see also, e.g., Benedict et al. 2002 Brandle et 
ar]|2012 ) . In the second pearls-on-a-strmg model, on the 



galaxy center, forming a smaller nuclear ring. This turns 
out entirely consistent with the observational result of 



other hand, star formation takes place preferentially at 
the contact points between a ring and dust lanes. This 
may happen because gas clouds with the largest densi 
ties are usually placed at the contact points due to orbit 



crowding (e.g., Ken ncy ct al. 1992; Rcynaud fc Downes 



1997[ [Kohno efar|il999; Hsich ct al. 20rT ). Since star 
clusters age as they orbit along the ring, this model nat- 
urally predicts a bipolar azimuthal age gradient of star 
clusters starting from the contact points (see also, e.g., 



Ryder ct al.'2001^'Allar d et al.|20 06','Maz zuca et al.|2008 
^kcr ct al. 2008; Ryde r et al.| |2010 van der Laan et al 
2013f . [Mazzuca et al.| (2008) found that - 50% of the' 
nuclear rings in their sample galaxies show azimuthal 
age gradients and that such galaxies have, on average, a 
larger value of the mean SFR than those without notice- 
able age gradients. 

Another interesting observational result concerns ra- 
dial locations of star clusters relative t o the nuclear rings 



20011 



In some galaxie s such as NGC 1512 (Maoz et al. 
and NGC 4314 ( [Benedict et al.|2002[ ), young star clusters 
are loc ated at larger radii th an the dense gas of nuclear 
rings. Martini et al. (|2003 ) also reported that out of 



123 barred galaxies in their sample, eight galaxies have 
strong nuclear spirals, all of which have star-forming re- 
gions outside the rings. By an alyzing multi-waveb and 
HST archive data of NGC 1672, [Jang fc Lee| ( [20131 ) re- 
cently identified hundreds of young and old star clusters 
with ages in the range ~ 1 — 10'^ Myr. They found that 
the clusters in the nuclear regions exhibit a systematic 
positive radial age gradient, such that older clusters tend 
to be located at larger galactocentric radii, farther away 
from the ring. Proposed mechanisms for the radial age 
gradient include the de crease in the ring size d ue to an- 
gular momentum loss (Regan & Teuben 20031) and mi- 
gration of clusters due to tid al interactions with the ring 
( [van de Ven fc C hang 2009). 

Numerical simulations have been a powerful tool to 
study formation and evolution of bar substructures such 
as dust lanes, nucle ar rings, and nuclear s pirals (e.g., 
Sanders fc Huntley|1 976; Athanassoula 1992; Finer et aL[ 



1995!'Englmai er fc Gerh ard 1997; Pat sis fc Ath anassoula' 
2000 ; Macicjewski et a l. 2002; Regan fc Teuben 2003i 
2004| [Ann fc Thakur[ [2005; Thakur ct al. 2009; KiS^ 



2012a|bp . In particular, ,Atlianassoula| ( ,1992, ) showed 



that dust lanes are shocks formed at the downstream side 
from the bar major axis. Dust lanes tend to be shorter 
and located close r to the bar major axis as the gas sound 
speed incre ases (Englmaier fc Gerhard 1997 Patsis fc 



Athanassoula„ 2000 ; Ki m et al. .. 2012a ~ 

Very recently, Kim ct aT ( |2012b[ hereafter Paper I) ran 
various models with differing bar strength and demon- 
strated that nuclear rings form not by resonant inter- 
actions of the gas with the bar potential, as was previ- 
ously thought, but instead by the centrifugal barrier that 
the inflowing gas with non- vanishing angular momentum 
cannot overcome. According to this idea, a more massive 
bar forms stronger dust-lane shocks which remove an- 
gular momentum more efficiently from the gas, so that 
the inflowing gas is able to move inward closer to the 



Comeron et al. (2010) that "stronger bars host smaller 
rings". Magnetic stress at the dust lanes takes away an- 
gular momentum additionally, leading to an even smaller 
ring compar ed to the unmagnetized counterpart (jKim fc] 



Stone 20121. Paper I also showed that nuclear spirals 
that form mside nuclear rings unwind with time due to 



the nonlinear effect ( Lee fc Goodman|1999 ), with an un- 
winding rate higher for a stronger bar. 'I'hus, the prob- 
ability of having more tightly wound spirals is larger for 
galaxies with a weaker bar, consistent w ith the observa- 
tional result of Peeples fc Martini (2006). 

While the numerical studies mentioned above are use- 
ful to understand gas dynamics in the central regions 
of barred galaxies, they are without self-gravity and/or 
prescriptions for star formation. There have been only a 
few numerical studies that considered star formation in 
nuclea r rings in a self-consistent way. [Heller fc Shlosman] 
( 1994 1 studied star formation in galactic disks that are 
unstable to bar formation. Using a smoothed particle 
hydrodynamics (SPH) combined with iV-body method, 
they found that star formation in barred galaxies occurs 
episodically, with a time scale of ~ 10 Myr, and that the 
associated SFR is well correlated with the mass accre- 
tion rate to the central b l ack h ole (BH) . These were con- 
flrmed by Knapen et al. ( 1995 ) who also found that tur- 
bulen ce driven by star torm ation tends to widen nuclear 



rings. Priedli fc Benz (1995') used another SPH-|-iV-body 



method to run various models with differing parameters, 
flnding that star formation in the nuclear regions first ex- 
perience s a burst phase before e ntering a quiescent phase 
(see also Martin fc Friedli 1997). Since these authors em- 
ployed a small number (~ lO'') of gas particles in their 
models, howe ver, they w e re un able to resolve the nuclear 
regions well. Kim et al. (2011) ran SPH simulations for 
star formation specific to the central mo l ecular zone in 
the Milky Way. While [Dobbs fc Pringle[ ( |2010[ ) studied 
cluster age distributions m spiral and barred galaxies, 
their results were based on SPH simulations that did not 
consider star formation and feedback. 

In this paper, we extend Paper I by including self- 
gravity and a prescription for star formation feedback. 
We focus on temporal and spatial distributions of star 
formation occurring in nuclear rings of strongly-barred 
galaxies. Unlike the previous SPH simulations with star 
formation, our models use a grid-based, cylindrical code 
with high spatial resolution in the central regions. We 
also allow for time delays between star formation and 
feedback, which is crucial to study age gradients of star 
clusters that form in nuclear rings. Our main objectives 
are to address important questions such as what controls 
the SFR in the nuclear rings and what are responsible 
for the presence (or absence) of the age gradients of star 
clusters in the rings, mentioned above. 

We take a simple galaxy model in which a self- 
gravitating gaseous disk with either uniform or expo- 
nential density distribution is placed under the influ- 
ence of a non-axisymmetric bar potential. We imple- 
ment a stochastic prescription for star formation that 
takes allowance for a threshold density as well as a star 
formation efficiency. Star formation feedback is treated 
only through direct momentum injections from super- 
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nova (SN) explosions occurring 10 Myr after star forma- 
tion events. By considering an isothermal equation of 
state, we do not consider gas cooling and heating, and 
radiative feedback, which may be import ant in regulat 



ing s t ar formation in disk galaxies (e.g., Ostrikcr et al. 



20T0l[Ostn 



Kim et al._ 2U11' rfhetiy 



cer & Stietty||2UlH Kim et al..,:^Ull, Stietty 
2012[ ). In our models, the bar potential is 



& Ostriker 

turned on slowly over time, which not only represents a 
situation where the bar forms and grows but also helps 
avoid abrupt gas responses. In each model, we measure 
the SFR in the ring and study its dependence on various 
quantities such as the gas mass in the ring, mass inflow 
rates to the central regions, bar growth time, etc. We also 
explore temporal and spatial variations of star-forming 
regions and their connection to the SFR. In addition, we 
study physical properties of star clusters and gas clouds 
in the rings and compare them with observational results 
available. 

We remark on a few important limitations of our mod- 
els from the outset. First of all, our gaseous disks are 
two-dimensional and razor-thin. This ignores potential 
dynamical consequences of vertical gas motions and re- 
lated mixing, which was shown important in ind ucing 



non-steady gas motions across spiral sho cks (e.g., Kim 



k Ostriker 2006 Kim et al. 2006 2010). SecondT"^ 



adopt an isothermal equation of state tor the gas, cor 
responding to the warm phase, and do not consider ra- 
diative cooling and heating required for produc t ion and 



transitions of multipha se gas (e.g., Field 



1965 Wolfire 



[et al. 2003^ McKce fc Ostriker] |2007l ). We also ignore 



the ettects of outflows, winds, and radiative feedback 
from young stars, which may be of crucial importance 
in setting up the equilibrium pressure in galactic planes, 
thereby reg ulating star formation in disk galaxie s (e.g., 
Ostriker et al. 2010; O striker fc Shetty 2011, Shetty fc| 



Ostriker 2012). Finally, we in the present work do not 
consider the effects of spiral arms that may supply gas to 
the bar regions. With these caveats, the numerical mod- 
els presented in this paper should be considered as a first 
step toward more realistic modeling of star formation in 
nuclear rings. 

This paper is organized as follows. In Section 2, we 
describe the numerical methods and parameters used for 
our time-dependent simulations. In Section 3, we present 
the temporal evolution of SFRs, and their dependence on 
the model parameters. The age gradients of star clusters 
in the azimuthal and radial directions as well as prop- 
erties of star clusters and dense clouds are discussed in 
Section 4. In Section 5, we summarize our main results 
and discuss their astronomical implications. 

2. MODEL AND METHOD 

To study star formation in nuclear rings of barred 
galaxies, we extend the numerical models studied in Pa- 
per I by including self-gravity, conversion of gas to stars, 
and feedback from star formation. In this section, we 
briefly summarize the current models and describe our 
handling of star formation and feedback. The reader is 
referred to Paper I for more detailed description of the 
numerical models. 

2.1. Galaxy Model 

We initially consider an infinitesimally-thin, rotating 
disk. The disk is assumed to be unmagnetized and 



Table 1 

Model Parameters 



Model So(M0pC-2) -TbarAo 



(1) 



(2) 



(3) 



(4) 



noSG 


20 




0. 


U05 


5 




0.75 


UIO 


10 




0.75 


U20 


20 




0.75 


U30 


30 




0.75 



M25 


20 


1 


0.25 


M50 


20 


1 


0.50 


FB05 


20 


0.5 


0.75 


FB20 


20 


2 


0.75 


FB40 


20 


4 


0.75 


E30 


30 


1 


0.75 


E50 


50 


1 


0.75 


ElOO 


100 


1 


0.75 



Note. — Gas surface density in Models 
with the prefix "E" initially have an exponen- 
tial distribution S = Sq exp (— r/3.5 kpc). All 
the other models have a uniform density dis- 
tribution E = Sq. 



isothermal with sound speed of Cs — 10 km s^^. The 
external gravitational potential responsible for the disk 
rotation consists of four components: a stellar disk, a 
stellar bulge, a non-axisymmetric stellar bar, and a cen- 
tral BH with mass Mbh = 4 x 10^ Mq. This gives rise to 
a rotation curve that is almost flat at Vc ~ 200 km s~^ 
in the bar region and its outside. The presence of the 
BH makes the rotation velocity rise as Vc oc {MBn/r)^^^ 
toward th e galaxy center . The bar potential is mod- 
eled by a Ferrers ( 1887 1 prolate spheroid with semi- 
major and minor axes of 5 kpc and 2 kpc. 



respec- 
tively. The bar is rigidly rotating with a pattern speed 
fib — 33 km s~^ kpc~ , which places the corotation res- 
onance radius at r = 6 kpc and the inner Lindblad res- 
onance (ILR) radius at r = 2.2 kpc. The corresponding 
orbital time is iorb — ^ir/il = 186 Myr. In our mod- 
els, the bar potential is turned on over the bar growth 
time scale Tbar, while the central density of the spheroidal 
component (bar plus bulge) is kept fixed. We vary Tbar 
to study situations where the bar grows at a different 
rate. The mass of the bar, when it is fully turned on, is 
set to 30% of the total mass of the spheroidal component 
within 10 kpc. All the models are run until 1 Gyr. 

As in Paper I, we integrate the basic equations of ideal 
hydrodynamics in a frame corotating with the bar. We 
use the CMHOG code in cylindrical polar coordinates 
(r, 0). CMHOG is third-order accurate in space and has 
very little numerical diffusion ( Piner et al.||1995" ). To re- 
solve the central region with high accuracy, we set up a 
logarithmically-spaced cylindrical grid over r = 0.05 kpc 
to 8 kpc. The number of zones in our models is 1024 in 
the radial direction and 632 in the azimuthal direction 
covering the half-plane from (p = — 7r/2 to tt/2. The cor- 
responding spatial resolution is 0.25 pc, 5 pc, and 40 pc 
at the inner boundary, at r = 1 kpc where most star 
formation takes places, and at the outer radial bound- 
ary, respectively. We adopt the outflow and continuous 
boundary conditions at the inner and outer radial bound- 
aries, respectively, while taking the periodic boundary 
conditions at (/) = ±7r/2. 

Our models consider conv ersion of gas to stars, as will 
be explained in Section [2?2] in detail. Since the total gas 
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mass transformed to stars is significant, it is important 
to evolve them under the combined gravitational poten- 
tial of the gas and stars. At each time step, we calcu- 
late the stellar surface density on the grid p oints via the 



triangular-shape d-cloud assignment scheme (Hockney & 
Eastwood 1988) from the distribution of stellar parti 



cies. We then solve the Poisson equation to obtain the 
gravitational potentia l of the total ( gas plus star) sur- 
fa ce density, using the Kalnajs] ( 1971 1 method presented 
Shetty fc Ostriker| ([2008) [J To allow for dilution of 



gravity due to tinite thickness H of the combined disk 
we take H/r = 0.1 as a softening parameter in the po- 
tential calculation. 

To explore the dependence of SFR upon the to- 
tal gas content and the way the gas is spatially dis- 
tributed, we initially consider gaseous disks with either 
uniform surface density Sq or an exponential distribution 
Eg exp(— r/i?d) with the scale length of Rd = 3.5 kpc. 
We also vary the fraction /mom of the radial momentum 
from SNe imparted to the disk in the in-plane direction 
relative to what would be the total radial momentum in 



a three-dimensional uniform medium (see Section 2.2). 
We run a total of 13 models that differ in Eq, Tbar, and 
/mom- Table [l] lists the model parameters. Column (1) 
lists each model. Models with the prefix "E" have an ex- 
ponential disk, while all the others have a uniform disk. 
Model noSG is a control model that does not include 
self-gravity and star formation. Column (2) lists Eq of 
the disk. Column (3) gives the bar growth time Tbar in 
units of iorb, while Column (4) gives /mom- We take 
Model U20 with Eq = 20 M© pc^^^ TbarAorb = 1, and 
/mom = 0.75 as our fiducial model. All the models ini- 
tially have a Toomre Q parameter greater than unity, so 
that they are gravitationally stable in the absence of a 
bar potential. However, nuclear rings that form near the 
center achieve large density, enough to undergo runaway 
collapse to form stars. 

2.2. Star Formation and Feedback 

To model star formation and ensuing feedback, we first 
identify high-density regions whose average surface den- 
sity (E) within a radius i?sF exceeds a critical density. 
The natural choice for the threshold density would be 



^th 



2GR 



•SF 



1160 M© pc" 



10 km s 



RsF 
10 pc 



(1) 
from the Jeans condition. While it is desirable to choose 

a small value for the sizes of star-forming regions, we 

take i?sF — 10 pc because of numerical resolution: a 

star- forming cloud at r ^ 1 kpc encloses typically '-^ 13 

grid points. 

Not all clouds with (E) > Eth immediately undergo 

gravitational collapse and star formation since we need 

to consider the star form ation efficiency as well as the 

computational time step (Kim et al. 2011). The SFR 



expected from a c loud with mass Mdoud 
from the Schmidt ( 1959 ) law, is 



^i?2p(S) 



SFR = eff 



M, 



cloud 



iff 



for (E) > E. 



th: 



(2) 



^ We ignore the gravity from the initial gas distribution in order 
to make the initial rotation curve the same with that in the non- 
self-gravitating counterpart. 



where eg is the star formation efficiency per free-fall time, 
iff, defined by 



*^ (32G(p)j 



1/2 



= 3.4 Myr 



(S) 



-1/2 



1160 Mq pc-2 

(3) 
assuming a disk scale height of 100 pc. We take eff — 
0.01, consistent with theoretical and observational esti- 



mates (e.g., Krumholz &: McKee[2005 Krumholz & Tan 
2007| ). 

The star formation probability of an eligible cloud with 
(E) > Eth in a time interval At i s then given by p = 
l-exp(-effAt/fff) w effAt/ts (e.g., [Hopkins et al.|20lT). 



For a given computational time step At, the probabUity 
p calculated in our models is typically ^ 10^^ — 10^^, 
much smaller than unity. In each time step, we thus 
generate a uniform random number Af € [0, 1), and turn 
on star formation only if A/" < p. When a cloud undergoes 
star formation, we create a particle with mass M*, and 
convert 90% of the cloud mass to the particle mass. The 
initial position and velocity of the particle are set equal 
to the density- weighted mean values of the parent cloud 
within Rsp . Each particle has a mass in the range M* ~ 
10^ — 10^ Mq, which is about ~ 1 — 10^ times larger 
than the masses of observed clusters in nuclear r ings (e.g 



200T1 [Benedict et al.|2002, see also |Portegiesl 



Maoz et al 

Zwart et al. 2010p . Therefore, a massive single particle in 
our models can be regarded as representing an unresolved 
group of star clusters rather than an individual cluster. 

We treat SN feedback using simple momentum input 
to the surrounding gaseous medium. We consider only 
Type II SN events since our models run only until 1 Gyr. 
Since we do not resolve individual stars in a cluster or 
their group, we assume that all SN explosions occur si- 
multaneously. Stars with mass between 8 M p and 40 Mq 
explode as Type II SNe ( Heger et al.|[2003 ), which com- 
prise about 7% of the cluster mass under the [Kroupa 
(2001 1 initial mass function. The mean mass of SN pro- 
gemtors is then ~ 14 M0 , indicating that the number of 
SNe exploding from a cluster (or their group) with mass 
M* is AsN — M*/(200 M©), with the total ejected mass 
-^^cjocta = 0.07Af* returning back to the ISM. Between 
star formation and SN explosions, we allow a time delay 
of 10 Myr, correspondi ng to the mean life time o f Type 
II SN progenitors (e.g., Lejeune fc Schaerer[[2001 ). Note 
that the orbital time or gas in nuclear rings is typically 
^ 25 Myr in our models, so that star clusters move by 
about 150° in the azimuthal angle from the formation 
sites before experiencing SN explosions. 

Each feedback, corresponding to A/sn simultane- 
ous SNe, injects mass and radial momentum in the 
form of an expanding shell. In the momentum- 

erg 
3 X 



conserving stage, a single SN with energy 10^^ 
would drive radial momentum Prad,3D = 
105eo^^(E/l M© pc-2)-i/4 Mq km s'^ to the surround- 
ing gas if the background medium is uniform and in three 
dimens ions, where ep is the SN energy in units of 10 ^^ erg 
(e.g., |Chevalier|1974[[Siuin|[r980, Cioffi et al.|1988[ ). The 



dependence ot i-Vad.sD on E is due to the tact that the 
shell expansion in the Sedov phase is slow when the back- 
ground surface density is large. Since our model disks are 
razor-thin by ignoring the vertical direction, the momen- 
tuin imparted to the gas in the simulation domain would 
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be smaller than -Piad,3D- Let /mom denote the fraction of 
the total radial momentum that goes into the in-pla ne di 
rectiqn. If the expansion is isotropic, /mom ~ 75% (Kim 
et al.|20lT ) , but /mom can be smaller in a vertically strat- 
ihed disk since it is easier for a shell to expand along the 
vertical direction. The total momentum of a shell from 
each feedback is thus set to 



P,h = 3xl0V, 



mom J^SN I ^ M^ p^_2 



-1/4 



Mp, kms"^ 

(4) 

In this paper, we take /mom = 0.75 as a standard value, 
but run some models with lower /mom to study the effect 
of /mom on the SFR (see Table m). 

As the initial radius of a shelX we take Rsh = 40 pc, 
corresponding to the shell size at the end of the Sedov 
phase when A/sn = lO'^ and the background density is 
E = 10^ Mq pc~^, the typical mean density of nuclear 
rings when star formation is active. When feedback oc- 
curs from a particle, we redistribute the mass and mo- 
mentum within a circular region with radius Rsh centered 
at the particle by taking their spatial averages. We then 
add the shell momentum density 



SVsh = 



Pn 
0, 



r > Rsh, 



(5) 



to the gas momentum density in the in-plane direction, 
and Aiojocta/(7ri?sii) to the gas surface density, while re- 
ducing the particle mass by Mojocta- In equation (|5|), R 
denotes the position vector relative to the particle loca- 
tion and Pmax = 2Psh/(7ri?g[j) is the momentum per unit 
area a,t r — Rsh- Note that Vsh{R) oc R^ ensures an 
initially divergcncc-frce condition at the feedback center 
(e.g., Kim et al.„2011J . 

3. STAR FORMATION IN NUCLEAR RINGS 

In this section, we first describe overall evolution of 
our numerical simulations. We then present the tempo- 
ral variations of SFRs and their dependence on the gas 
mass, the bar growth time, and /mom, as well as the re- 
lationship between the SFR surface density and the gas 
surface density. The properties of star clusters and gas 
clouds in the rings are analyzed in the next section. 

3.1. Overall Gas Evolution 

We begin by presenting evolution of our standard 
model U20 that has a uniform density So = 20 Mq pc~^ 
initially. Evolution of the other models is qualitatively 
similar, although more massive disks show more active 
star formation. In all models, star formation occurs 
mostly in the nuclear rings jj Figure [l] plots snapshots 
of the gas distributions as well as the positions of star 
clusters in Model U20 at a few selected epochs. The 
left panels show the gas density in logarithmic scales 
in the 5 kpc regions, with the solid ovals indicating 
the outermost xi-orbit that cuts the x- and y-axes at 
Xc = 3.6 kpc and yc = 4.7 kpc, respectively. The middle 
and right panels zoom in the central 2 kpc regions. In the 
right panels, small dots indicate star clusters older than 

^ In Model U30, the bar-end regions, often called ansae, form 
stars as well, although the associated SFR is less than 3% compared 
to the ring star formation. 



10 Myr, while asterisks correspond to young clusters with 
age < 10 Myr: the upper right colorbar represents their 
ages. In all panels, the bar is oriented vertically along 
the y-axis and remains stationary. The gas inside the 
corotation resonance is rotating in the counterclockwise 
direction. 

Early evolution of Model U20 before t = 0.1 Gyr 
is not much different from non-s elf-gravitating models 
presented in 
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( |2012a[ ). The imposed non- 
axisymmetric bar potential perturbs gas orbits to cre- 
ate overdense ridges at the far downstream side from the 
bar major axis. Only the region inside the outermost xi- 
orbit responds strongly to the bar potential. Before the 
bar is fully turned on (i.e., t < Tbar), no closed xi- and 
a;2-orbits exist since the external gravitational potential 
varies with time. During this time, the gas streamlines 
are highly transient as they try to adjust to the time- 
varying potential. The overdense ridges grow with time 
and move slowly toward the bar major axis, develop- 
ing into dust-lane shocks. The gas passing through the 
shocks loses angular momentum and flows radially in- 
ward along the dust lanes. The radial velocity of the 
inflowing gas is so large that it is not halted at the ILR 
(Paper I). It gradually rotates faster due to the Goriolis 
force, and eventually forms a nuclear ring after hitting 
the dust lane at the opposite side. Produced by super- 
sonic collisions of two gas streams, the contact points 
between the ring and the dust lanes have the largest den- 
sity in the ring, and can thus be preferred sites of star 
formation. 

Over time, the nuclear ring shrinks in size and the 
contact points rotate in the azimuthal directionP] To 
illustrate this more clearly. Figure [2l plots the temporal 
and radial variations of the azimuthally averaged surface 
density in the central regions of Model U20 and its non- 
self-gravitating counterpart. Model noSG. The horizontal 
dashed line marks the location of the ILR. The ring is 
beginning to form at t ~ 0.1 Gyr at the galactocentric 
radius of r ~ 1 kpc, well inside the ILR. At early time 
when the bar is growing, the ring is not in an equilibrium 
position and its shape is quite different from a:;2-orbits: 
the major axis of the ring is inclined to the x-axis by 
-30° at t = 0.15 Gyr (Fig. [ik). As the ring material 
continuously interacts with thebar potential, the contact 
points rotate in the counterclockwise direction. At i ~ 
0.2 Gyr, the ring settles on one of the j:2-orbits, and the 
contact points are located close the bar minor axis. 

At the same time, the ring becomes smaller in size due 
to the addition of low-angular momentum gas from out- 
side as well as by collisions of the ring material whose or- 
bits are perturbed by thermal pressure (Paper I). When 
self-gravity is absent, the decreasing rate of the ring ra- 
dius is dlnrNRM 1-5 Gyr"^ untfl t < 0.8 Gyr. Af- 
ter this time, the ring is so small that strong centrifu- 
gal force inhibits further decay of the ring. In Model 
U20 with self-gravity included, on the other hand, strong 
self-gravitational potential of the ring makes the gas or- 

^ As pointed out by the referee, Figure 111 shows that the ring 
starts out fairly elongated and titled relative to the bar and be- 
comes rounder with time. Note that the shape of the ring in 
Model U20 at t = 0.15 Gyr is remarkably similar to t hat of a 
highly-elo ngated nuclear ring in ESO 565-11 observed by |Buta et| 
|al.| ( |1999| l, although the latter is ^ 2-3 times bigger than the tor- 
mer. 
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Figure 1. Snapshots of logarithm of gas density (color scale) as well as the locations of star clusters in Model U20 at i = 
0.15,0.2,0.5,0.8 Gyr. The left panels show the 5 kpc regions, while the middle and right panels zoom in the central 2 kpc regions. 
The upper left colorbar labels log(E/So). The solid ovals in the left panels draw the outermost xi-orbit that cuts the x- and y— axes at 
Xc = 3.6 kpc and j/c = 4.7 kpc. The dotted curve in (d) is an xi-orbit with Xc = 1.7 kpc and pc = 4.4 kpc that traces the inner ring. Small 
dots in the right panels denote clusters older than 10 Myr, while asterisks represent clusters younger than 10 Myr, with the upper right 
colorbar displaying their ages. 
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Figure 2. Temporal and radial variations of the azimuthally av- 
eraged surface density in logarithmic scale for (a) Model U20 and 
(b) Model noSG. The horizontal dashed line in each panel marks 
the location of the inner Lindblad resonance. The colorbars label 
log(E/So). The decreasing rate of the ring size is smaller when 
self-gravity and star formation are included. 

bits relatively intact, reducing the ring decay rate to 
d\nrjssR/dt ~ —0.4 Gyr~^. The decrease of the ring size 
in turn causes the contact points to move radially inward 
and rotate further in the counterclockwise direction (but 
not more than ~ 30°). 

As Figure ll] shows, the bar region (i.e., inside the 
outermost-xi orbit) becomes evacuated quite rapidly. 
This is because the bar potential is efficient in remov- 
ing angular momentum from the gas only in the bar re- 
gion, while its influence on the gas orbits outside the 
outermost-xi orbit is not significant. By t ~ 0.3 Gyr, 
most of the gas in the bar region is transferred to the 
ring. The amount of the gas added to the bar region 
from outside is much smaller than that lost to the ring, 
which causes the mass inflow rate to the ring to decrease 
dramatically with time (see Section 3.2). 

Figure [TJi shows that at late time there is a substan- 
tial amount of gas trapped in around an xi-orbit with 
Xc = 1.7 kpc and i/c — 4.4 kpc, shown as a dotted line. 
This elongated gaseous structure circumscribing the dust 
lanes and nuclear ring is called the inner ring (e.g., Buta 
[1986; Reg an et al.|2002 ): we term the corresponding Xi- 
orbit the inner-ring xi-orbits. The formation of the inner 
ring in our model is as follows. As mentioned before, the 
dust-lane shocks form first at far downstream and moves 
toward the bar major axis as the bar potential grows. 
During this time, much of the gas in the bar region in- 
falls to the nuclear region. Near the time when the bar 
attains its full strength, the dust lanes find their equi- 
librium positions on an a;i-orbit, still at the leading side 
from the bar major axis. The inner ring starts to form 
at this time, by gathering the residual material located 
between the outermost and inner-ring xi-orbits. Some 
gases located outside the outermost a;i-orbit experience 
collisions near the bar ends where a:i-orbits crowd, and 
are then able to lower their orbits to the inner-ring xi- 
orbit, increasing the inner-ring mass. In Model U20, the 
gas added to the inner ring from outside of the outer- 



most xi-orbit is about 70% of the total inner-ring mass 
at t — 1 Gyr. Most of the gas inside the inner-ring 
a; 1 -orbit had already transited to the nuclear ring by ex- 
periencing the dust-lane shocks before the bar was fully 
turned on. 

3.2. Star Formation Rate 

Figure [3] plots the time evolution of the SFR, the mass 
inflow rate to the nuclear ring Mnr = L TjVrrdcj) (mea- 
sured at r — 1.5 kpc), and the total gas mass inside 
the nuclear ring Mnr in Model U20. Here Vr denotes 
the radial velocity of the gas. In plotting these profiles, 
we take a boxcar average, with a window of 20 Myr. In 
this model, the bar potential grows over the time scale of 
■''bar = 0.19 Gyr, and the first star formation takes place 
at the contact points at t = 0.12 Gyr. As the bar grows 
further, the dust-lane shocks become stronger, increasing 
the amount of the infalling gas to the ring. Mj^n attains 
a peak value ~ 8 Mq yr~^ at t = 0.15 Gyr, which coin- 
cides with the time of highest density of the dust lanes 
(see Fig. 5 of Paper I) . The decay of Mnr after the peak 
is caused by the fact that only the gas inside the out- 
ermost xi-orbit responds strongly to the bar potential, 
while the outer region is not much affected. Similarly, 
the SFR exhibits a strong burst with a maximum value 
~ 8 M0 yr~^, which occurs ~ 30 Myr after the peak of 
Mnr. The associated SN feedback produces many holes 
in the gas distribution, driving a huge amount of kinetic 
energy to the surrounding medium. Note that the nu- 
clear ring, albeit somewhat patchy, is overall well main- 
tained despite energetic momentum injections (Fig. Hp). 

When the mass inflow rate to the ring is very large, star 
formation occurring at the contact points alone is unable 
to consume the whole inflowing gas. As we will show in 
Section [4.1[ the maximum gas consumption rate afforded 
to the contact points is estimated to be about 1 M0 yr~^ 
for the parameters we adopt. Any surplus inflowing gas 
passes by the contact points and is subsequently added 
to the ring that is clumpy. Some overdense regions in 
the ring are soon able to achieve the mean density larger 
than the critical value and undergo star formation, in- 
creasing the SFR rapidly. Since the star-forming regions 
are randomly distributed in the ring, there is no obvious 
azimuthal age gradient of star clusters in this high-SFR 
phase. 

Particles spawned from star formation orbit about the 
galaxy center under the total gravity, but they do not feel 
gas pressure that is quite strong in the nuclear ring. In 
addition, the gaseous ring becomes smaller in size with 
time. Thus, the orbits of star particles increasingly de- 
viate from the gaseous orbits over time. When stars in 
clusters explode as SNe, they are not always located in 
the ring. A majority of star clusters are still in the ring, 
while there are some clusters (~ 10%) located exterior to 
the ring. SNe occurring in the ring have relatively small 
radial velocities due to a large background density, and 
the shell expansion is limited by the surrounding dense 
gas. On the other hand, SNe exploding outside the ring 
can have very large expansion velocities enough to send 
the neighboring gas out to the bar-end regions. The ex- 
pelled gas sweeps up the gas on its way to the bar-end 
regions, passes through the dust-lane shocks again, and 
falls back radially inward along the dust lane. This in- 
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Figure 3. Temporal variations of (a) tlie SFR, (b) tlie mass inflow rate Mjmr to tlie nuclear ring, and (c) the total gas mass Mnr in the 
nuclear ring of Model U20. The ordinate is in linear scale in the left panels, while it is in logarithmic scale in the right panels. In (a) and 
(b), the horizontal dotted lines indicate the reference rate of 1 TsAq yr~^. 



creases Mnr temporarily during t — 0.2-0.35 Gyr, with 
the associated short bursts of star formation a,t t = 0.27 
and 0.32 Gyr in Model U20. 

At i = 0.3 Gyr, the bar region is almost emptied, ex- 
cept for the inner ring, as most of the gas is already lost 
to the ring. Other than intermittent infalls of the ex- 
pelled and swept-up gas from SNe, the mass inflow rate 
from the inner ring to the nuclear ring and the related 
SFR become fairly small. In addition, Mnr is reduced 
to below 2 X 10^ Mq in Model U20 since star formation 
consumes the gas in the ring, which also decreases the 
SFR (Fig. [3]). As most of the gas in the bar region inside 
the outermost a;i-orbit is almost lost to star formation, 
the galaxy evolves into a quasi-steady state where star 
formation is limited to small regions near the contact 
points (Fig.[l}i). 

3.3. Parametric Dependence of SFR 

Figure [4] compares the SFRs from (left) uniform-disk 
and (right) exponential-disk models with different Eq. 
In all models, the SFR displays a primary burst fol- 
lowed by a few secondary bursts, with time intervals 



of ~ 50 — 80 Myr, before becoming reduced to below 
1 Mq yr~^. The primary burst is associated with the 
rapid gas infall due to angular momentum loss at the 
dust-lane shocks, while the secondary bursts are caused 
by the re-infall of the ejected gas via SN feedback out 
to the bar region. Models with larger Eq start to form 
stars earlier and have a larger value of the maximum 
star formation rate, SFR^ax- The duration of active 
star formation, Aigp, defined by the time span when 
SFR > SFRniax/2, is also larger for models with larger 
Eo, since the gas available for star formation is corre- 
spondingly larger. Models U20 and E50 initially have 
a similar gas mass inside the outermost a;i-orbit, but 
Model E50 has larger SFRmax since the gas is more cen- 
trally concentrated and thus infalls more readily to the 
nuclear ring. Columns (2) and (3) of Table [2] list SFR^ax 
and AigF for all models. The phase of active star forma- 
tion lasts only for ^ 0.5Tbar in all models. 

Figure [s] shows how (left) the bar growth time and 
(right) the amount of the momentum injection affect the 
temporal behavior of the SFR. In models where the bar 
grows more rapidly, dust-lane shocks form earlier and 
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Figure 5. Temporal variations of the SFR for (a) models with differing the bar growth time T^ar and (b) models with different momentum 
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Table 2 

Simulation Outcomes 



Model 


SFRinax 


AtsF 


Ml°^ 


A/,1 


Mnr 


P 


r 




(M0 yr-i) 


(Myr) 


(1O«M0) 


(lO^Mo) 


(1O«M0) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(S) 


noSG 


_ 


_ 


_ 


10.66 


8.03 


_ 


_ 


U05 


0.53 


- 


1.4 


2.62 


0.87 


13.7 


3.1 


UIO 


3.55 


42 


3.3 


5.33 


1.39 


14.4 


2.6 


U20 


7.46 


83 


9.0 


10.66 


1.42 


12.5 


2.2 


U30 


11.8 


83 


13.3 


15.99 


2.36 


12.6 


2.2 


M25 


7.77 


87 


8.2 


10.66 


1.08 


9.6 


2.3 


M50 


7.78 


77 


8.7 


10.66 


1.55 


10.8 


2.2 


FB05 


9.83 


68 


8.9 


10.66 


1.62 


12.3 


2.2 


FB20 


4.04 


163 


8.7 


10.66 


1.73 


9.7 


2.2 


FB40 


1.94 


392 


7.2 


10.66 


2.75 


7.4 


3.0 


E30 


5.10 


74 


6.1 


7.55 


1.55 


12.4 


2.2 


E50 


10.4 


77 


10.1 


12.59 


1.41 


6.3 


2.2 


ElOO 


21.9 


84 


23.0 


25.18 


1.78 


6.3 


2.0 



Note. — SFRjnax and Atgp denote the peak rate and the duration of active 
star formation with SFR > SFRinax/2, respectively; M*°* is the total mass in 
stars at t = 1 Gyr; Mx^ is the total gas mass inside the outermost xi-orbit at t = 
0; Mnr is the mass of the nuclear ring at t = 1 Gyr; /3 = rflog(t/yr)/d(r/kpc) 
is the radial age gradient of clusters; F = — dlog A'^/dlogM* is the power-law 
slope of the cluster mass functions. 
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Figure 6. Dependence of the SFR on (a) the mass inflow rate A/nr to the ring and (b) the total mass A/nr in the ring. The dashed line 
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correlated with A'/nr when SFR < 1 Mq yr"-"^. 
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Figure 7. (a) Dependence of the SFR surface density Sspa 
on tile mean surface density E^i of dense ciouds for Modei U20. 
Tlie colorbar indicates tiie epocli of star formation for eacli sym- 
bol, and tlie curved arrow denotes tlie mean evolutionary track 
in the SgFR — Sd plane, (b) Ssfr — Scl relationship from our 
models compared to the observed Kennicutt-Schmidt law. Squares 
and pluses repr esent norrnal an d circumnuclear starburst galaxies 
adopted from , Kennicutt| l |1998| l , respectively, while circl es are for 
spatiallv- resolved star-torming regions in NGC 1097 from |Hsieh et| 

p^ponl i. 

initiate stronger gas inflows. This causes star formation 
in models with smaller Tbar to occur at a higher rate and 
for a shorter period of time (see Table [2]). The peak SFR 
is attained approximately at i ^ (0.8 — l)rbar- This sug- 
gests that galaxies in which the bar forms more slowly are 
likely to have star formation less active instantaneously 
but extended for a longer period of time. FigurelSb shows 
that the SFR computed in our models is largely insensi- 
tive to /mom, although smaller /mom makes the secondary 
bursts less active. 

To directly address what controls star formation in nu- 
clear rings, we plot in Figure|6]the dependence of the SFR 
(left) on the mass inflow rate to the ring and (right) 
on the gas mass in the ring for Models UIO, U20, and 
U30. Color indicates the star formation epoch of each 
symbol. Even though there are large scatters especially 
when the SFR is low, the SFR is almost equal to Mnr 
over two orders of magnitude variations in Mnr. The 
scatters in the SFR — A/nr relation are due to the fact 
that star formation is stochastic in our models and that 
it takes the gas some finite time (~ 10 — 30 Myr) to travel 

from r = 1.5 kpc (where A/nr is measured) to the nu- 
clear ring. On the other hand, the SFR does not show 
a good correlation with A/nr. While SFR ex M^-^ for 
SFR /t 1 M0 yr^^, it is almost independent of A/nr for 
SFR ;$, 1 Mq yr^^. Note that the change in A/nr is less 
than a factor of 5 in Figure [SJd, while the SFR varies by 
more than two orders of magnitude. This suggests that 
it is the mass inflow rate to the ring, rather than the 
ring mass, that determines the SFR in the nuclear ring. 
Conversely, the SFR can be a good measure of the mass 



inflow rate driven by the bar potential. 

Column (4) of Table [2] gives the total mass in stars 
^^tot formed until the end of the run for each model. 
Columns (5) and (6) list the total gas mass M^-^ inside 
the outermost xi-orbit in the initial disk and the mass 
of the nuclear ring A/nr at t = 1 Gyr, respectively. Note 
that Mx-^ is approximately the maximum gas mass avail- 
able for star formation in the ring. We flnd that the 
relation Af*°* = M^^ - A//nr, with AfNR = 2 x 10® M© 
fixed, explains the numerical results fairly well, indicat- 
ing that most of the gas inside the outermost a;i-orbit 
flows inward to form stars, with some residual gas re- 
maining in the nuclear ring. Compared to Model U20 
with Tbar/iorb = 1, Modcl FB40 with Tbar/^orb — 4 has 
Af *°' about 20% smaller, since the gas in the bar region 
is still flowing in to the nuclear ring at the end of the 
run. 

As will be discussed in more detail in Section [5?2| the 
overall temporal trend of the SFR (that is, rapid de- 
cline after a primary burst except for a few secondary 
bursts) found in our numerical models is largel y similar 
to the nunierical results of previou s studies (e.g., HcUcr &| 
' Shlosman|1994[|Knapen et al.ll995 ; Friedli fc Bcnz 1995,), 



but appe ars m consistent w i th ob servations ot |Allard et| 



all ( 2006 ) and Sarzi et al. ( 2007 1 who found that star- 
lormmg nuclear rings live long, with multiple episodes of 
starburst activities. Since the ring SFR is controlled by 
the mass inflows rate to the rings, this implies that rings 
in real galaxies should be continually supplied with fresh 
gas from outside for quite a long period of time. Candi- 
date mechanisms for additional gas inflows, over a time 
scale much longer than the bar growth time, include spi- 
ral arms and cosmic gas infalls, which are not included 
in this paper. 

3.4. Star Formation Law 

To explore the dependence of the local SFR on the 
local gas surface density, we define clouds as regions 
in the simulation domain whose density is larger than 
300 M0 pc~^. This density roughly corresponds to 
the mean density of boundaries o f gra vitationally bound 
clouds in our models (see Section 4.4). While this choice 
of the minimum density for clouds is somewhat arbi- 
trary, these clouds may represent giant molecular clouds 
and their complexes includ ing hydrogen envel opes (e.g., 
Wniiams et aL|[2000{ |McKee fc Ostrike r 2007). 

At each time, we calculate the mean surface density 
Eel of, and the total area Ac\ occupied by, the clouds 
distributed along the ring. The mean SFR surface den- 
sity is then given by Esfr = SFR/Ad. Figure [7^ plots 
the resulting Esfr as a function of Sd for Model U20. 
Color represents the time when each point is measured, 
while the curved arrow indicates the mean evolution- 
ary direction in the Esfr~Sci plane. When the first 
star formation takes place {t — 0.12 Gyr), the ring has 
Eel - 560 Mq pc-2 and Ssfr ~ 0.5 M© yr"! kpc"^. 
The radial gas inflow along the dust lanes increases Esfr 
rapidly until it achieves a peak value at t = 0.18 Gyr. 
The corresponding increase of Ed is smaller since star 
formation reduces the gas content in the ring. After the 
peak, EgpR decreases with decreasing AfNR, but it has a 
larger value, by a factor of ^ 2 — 4 on average, than that 
at the same Ed before the peak. This is because active 
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Figure 8. Histograms of the star clusters that formed in each selected 
as a function of the azimuthal angle where they form. The arrow at 
point. 

SN feedback stirs the ring material vigorously, tending to 
increase density contrast between clumps and the back- 
ground material. After the Esfr peak, therefore, the 
total area covered by the gas with S > 300 Mq pc~^ 
becomes smaller than before, resulting in larger Ssfr- 

Figure [7p plots the Ssfr-Sci relationship measured 
at every 0.1 Gyr for (diamonds) the uniform-disk mod- 
els and (asterisks) the exponential-disk models. Models 
with larger M^-^ tend to have larger Esfr and larger 
Eel. Note that our numerical results are overall consis- 
tent with the Kennicutt-Schmidt law for normal galaxies 
(squares) and circumnucl e ar sta rburst galaxies (pluses) 
adopted from Kennicutt (1998), and not much differ- 
ent from the observed Esfr-Eci relation for spatially- 
resolved star-forming regio ns (circles) in the nu clear ring 

" There 
orders of 



(|2011|)F1 
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of NGC 1097 taken from [Hsieh et al. 
are large scatters in Esfr, amountmg to 
magnitude, both in observational and simulation results. 
For star formation in nuclear rings of barred galaxies, 
the gas surface density probably sets the mean value of 
Esfr, as the Kennicutt-Schmidt law implies, while the 
scatters in Esfr are likely due to the temporal variations 
of the mass inflow rate to the ring. 

4. PROPERTIES OF STAR CLUSTERS AND GAS CLOUDS 
4.1. Azimuthal Age Gradient 



time bin (with bin width of 0.15 Gyr) for all uniform-density models 
the bottom of each panel indicates the mean position of a contact 

As mentioned in Introduction, observations indicate 
that some galaxies have well-defined azimu thal age gra- 
dients of star cluste rs in nuc lear rings (e.g., 'R yder et al 
pOOTj lAUard et al.ir 2006; Bo ker et al. ;2008; R yder et al 
12010 ^ [van der L aan ct al 201 3]), while othe rs do not (e.g., 
|Benedict et al.||2002j Rrandle et al. 2012). Our simula- 
tions show that the presence or absence oi the azimuthal 
age gradient is decided by the SFR in the ring (or, more 
fundamentally, on Mnr), independent of Eg and the ini- 
tial gas distribution. 

Figure Is] plots for the uniform-density models the his- 
tograms of star clusters formed in each selected time 
bin, with bin width of 0.15 Gyr, as a function of the 
angular position where they form. The arrow at the 
bottom of each panel marks the location of a contact 
point that is moving in the positive azimuthal direction 
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- Mo 



Model U05 with 
1 (Fig. 11 



yr 



* In plotting the EgpR-S relation, [Hsieh et al.| ([201l| used the 
maximum density of a cloud, instead ot tne mean density, for E. 



with time, as described in Section 
Eo = 5 M© yr-i has SFR < _ 

and star-forming regions in this model are almost a 
ways localized to the contact point. In Models U20 and 
U30, on the other hand, star-forming regions are widely 
distributed along the azimuthal direction at early time 
(i = 0.15 - 0.45 Gyr) when SFR > 1 M© yr^^, while 
they are preferentially found near the contact points at 
late time {t > 0.45 Gyr) when SFR < 1 Mq yr'^. 

Star clusters age as they orbit along the ring and emit 
copious UV radiations during about ^ 10 Myr after 
birth. If star-forming regions arc localized to the con- 
tact points, therefore, clusters would appear as "pearls 
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Figure 9. Positions and ages of young star clusters at selected epoches in Models U05, U20, and U30, overlaid on the gas density 
distribution in linear scale. The left panels show a clear azimuthal age gradient, while there is no age gradient in the right panels. Color 
indicates the cluster ages. 

ure [9] for Model U30 t = 0.33 Gyr and Model U20 at 
t = 0.18 Gyr. 

Why does then the SFR (or Mnr) matter for the az- 
imuthal distributions of star-forming regions? The an- 
swer lies at the fact there is a limit on the rate of gas 
consumption at the contact points that occupy very small 
areas in the ring. When Mnr is sufficiently small, most 
of the inflowing gas to the ring can be converted into 
stars at the contact points, and the resulting SFR is cor- 
respondingly small. When Mnr is very large, on the 
other hand, the contact points cannot transform all the 
inflowing gas to stars instantaneously. The excess inflow- 
ing gas overflows the contact points and is transferred to 
other regions in the ring. The ring becomes denser not 



on a string" (e.g., Boker et al. 2008), with an age gra- 
dient along the rotational direction of the nuclear ring. 
This is illustrated in the left panels of Figure |9] which plot 
the spatial locations of star clusters with color indicating 
their ages (< 10 Myr), overlaid on the density distribu- 
tion in linear scale, in Model UIO at i = 0.35 Gyr and 
Model U20 at i = 0.51 Gyr. There is clearly a positive bi- 
polar age gradient starting from the contact points that 
are located at ~ 30° and 210°. On the other hand, 
when the star-forming regions are randomly distributed 
throughout t he ring, as in the "popcorn" model of Boker 
et al. ( 2008[ ), star clusters with different ages would be 
mixed. In this case, there is no apparent age gradient 
along the ring, as exemplifled in the right panels of Fig- 
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only by the addition of the overflowing gas but also by 
its own self-gravity. Some clumps in the ring achieve sur- 
face density above the threshold value, and start to form 
stars. We find that the rings have the Toomre stabil- 
ity parameter as low as ^ 0.5 when the SFR is near its 
peak, suggesting that gravitational instability promotes 
star formation. 

Observations show that galaxies with no azimuthal age 
gradie nt have, on average , slightly larger SFRs in the 
rings (Ma zzuca et al.|[2008 ), although the spread in SFR 
for individual galaxies is too large to make this conclu- 
sive. Note, however, that the critical SFR determining 
the presence or absence of the azimuthal age gradient de- 
pends on many parameters. The maximum SFR at the 
contact points can be estimated as follows. Let Ar and 
A0 denote the radial thickness and the azimuthal extent 
of a contact point, respectively. Then, the maximum 
SFR expected from two contact points is simply 



Af*,CP = 2effi;cprNRArA0/tff, 



(6) 



where Ecp is the surface density of the contact 
points. The mean density of star-forming clouds is 
- 4000 M© pc~^ which we take for Ecp. For eg = 0.01, 
Ar = 50 pc, A0 = 30°, and rNR = 1 kpc typical in our 
models, equation M yields Af*,cp ^ 1 M© yr~^, consis- 
tent with our numerical results. Note that the specific 
value of M*,cp depends on the parameters we adopt. In 

particular, M*,cp oc Cgr^j^/i?gp if the ring width is pro- 
portional to the ring radius, suggesting that galaxies with 
a weak bar (to have a smaller ring) and strong turbulence 

would have large Af^^cp- 

4.2. Radial Age Gradient 

Although the presence of an azimuthal age gradient 
depends on the SFR, we find that star clusters always 
display a radial age gradient. Figure 10 plots the spatial 
distributions of (left) the formation locations and (right) 
the present positions of star clusters on the x-y plane 
at i = 1 Gyr in Model U20. Each cluster is colored 
according to its formation time. The open circles denote 
the clusters that have passed the central region with r < 
0.3 kpc at least once during their orbital motions: such 
clusters would have been destroyed at least partially by 
strong tidal fields near the galaxy center if their internal 
evolution such as core collapse, evaporation, disruption, 
etc. had been considered. On the other hand, the filled 
circles are for clusters that have never approached the 
central region, and thus are most likely to survive the 
galactic tide. The dashed lines draw the ring at i = 
1 Gyr. Star-forming regions at late time are concentrated 
on the contact points, while they are well distributed at 
early time. Note that clusters that form early before 
the nuclear ring settles on an X2-orbit have initial kick 
velocities quite different from those on xi- or a:2-orbits at 
their formation locations. Although the ring soon takes 
on the a;2-orbit, these clusters move on eccentric orbits 
and wander around the nuclear region. Figure [T0| 3 shows 
that young clusters are preferentially found near the ring, 
while old clusters are located away from it, indicative of 
a positive radial gradient of their ages. 



(circles) at i — 1 Gyr as well as their formation loca- 
tions (plus symbols) for Models UIO and U20. Again, 
the open circles denote the clusters that have passed by 
the galaxy center, while the filled circles are for those 
that have not. Note that the age distributions of the 
present-day and formation-epoch locations of the clus- 
ters are not much different from each other, although 
the former shows a large spatial dispersion. The dis- 
persion is larger for older clusters. To quantify the ra- 
dial age gradient, we bin the clusters according to their 
ages, with a bin size of Alog(t/yr) = 0.2, and calcu- 
late the mean age and position in each bin. Our best 
fits of the ages to the formation-epoch positions are 
/3 = <ilog(t/yr)/d(r/kpc) - 14.4 and 12.5 for Models 
UIO and U20, respectively. Column (7) of Table [2] gives 
f3 for all models. This radial age gradient results pri- 
marily from the decrease in the ring size with time, such 
that old clusters formed at larger galactocentric radii. 
Clusters diffuse out radially through gravitational inter- 
actions themselves and also with dense clouds in the ring, 
without much effect on the radial age gradient. 

4.3. Cluster Mass Functions 

Figure [12] plots the mass functions of all the clusters 
that have formed in each of the (left) uniform-disk and 
(right) exponential-disk models until t — I Gyr. The 
upper panels are for the clusters formed while star for- 
mation is very active with SFR > 1 M© yr~^, whereas 
those produced when SFR < 1 Mq yr~^ are presented 
in the lower panels. In general, the mass distribution of 
clusters is described roughly by a power law, with its in- 
dex depending on the SFR and Afa;i . When the SFR is 
larger than 1 Mq yr~^, the slope of the mass function is 
F = —d log N/d log Af* ~ 2-3, with a smaller value corre- 
sponding to larger M^^- When SFR < 1 Mq yr"-*^, clus- 
ters have a much steeper mass distribution with F ;J, 3. 
Column (8) of Table [2] gives F. 

The dependence ofthe power-law index on the SFR 
and Afa;j can be understood as follows. When the SFR 
is large, there are numerous dense regions distributed 
throughout the ring. Such regions grow by accreting 
the surrounding material. Since star formation occurs 
in a stochastic manner in our models, some dense clouds 
have a chance to grow as massive as, or even larger than 
^ 10^ Mq, leading to a relatively shallow mass function. 
On the other hand, when the SFR is small, there are 
not many dense regions. Since the growth of density is 
quite slow in these models, they form stars at densities 
slightly above the critical value. In this case, most clus- 
ters have mass around ~ 10^ — 10^ M©, with a steep 
mass distribution. 



4.4. Giant Glouds 

Finally, we present the properties of giant clouds lo- 
cated in nuclear rings. High-resolution radio observa- 
tions show that nuclear rings consist of giant molecular 
associations at scale of ~ 0.2-0.3 kpc in which most star 
formation takes place. They typically have niasses of 



10 M f7^ and are gravitationally bound (e.g., Hsieh et 



l20i 



al. 2011) 



To show this more clearly. Figure 11 plots the age 
of clusters as functions of their present radial positions 



3entify giant clouds in our models, we utilize acore- 
finding technique developed by Gong & Ostriker (2011 ). 
This method makes use of the gravitational potential of 
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Figure 10. Spatial distributions of (a) the formation locations and (b) the present positions of star clusters in Model U20 at t = 1 Gyr. 
The dashed lines draw the ring at this time. Open circles denote the clusters that have passed the central region with r < 0.3 kpc during 
their orbits, while filled circles represent those that have not. Colorbar indicates the formation epoch of the clusters in unit of Gyr. At late 
time, stars form preferentially near the contact points. Gravitational interactions lead to diffusion of the clusters. 
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Figure 11. Ages of star clusters as functions of their current 
radial locations (circles) at t = 1 Gyr and their formation positions 
(pluses) for Models UIO and U20. Open circles are those that have 
passed by the galaxy center at a very close distance during orbital 
motions, while the filled circles are for those that have not. The 
dashed lines are the fits, with slopes of /3 = dlog(t/yr)/d(r/kpc) = 
14.4 and 12.5, for the initial cluster positions, for Models UIO and 
U20, respectively. 



the gas, and thus allows smoother cloud boundari es than 
the rriethod s based on isodensity surfaces (see, e.g., Smith 
et al. 2009). At a given time, we search for all the local 
mmirna ot the gravitational potential and find the largest 
closed potential contour encompassing one and only one 
potential minimum. We then define the potential min- 
imum and outermost contour as the center and bound- 
ary of a cloud, respectively. If the distance between two 
neighboring minima is less than 0.1 kpc, we combine 
them. Figure [13] plots, for example, giant clouds iden- 
tified by this technique for Model U20 at i = 0.36 Gyr. 
The left panel shows the gas surface density in linear 
scale, while the right panel displays boundaries of giant 
clouds as contours overlaid over the gravitational poten- 
tial of the gas. A total of 14 giant clouds are identified. 
The mean values of their masses M, radii i?, and one- 
dimensional velocity dispersions a are 10^ Mq, 100 pc, 
and 20 km s~^, respectively, corresponding to super- 
sonic internal motions. The average density of the cloud 
boundaries is found to be ~ 300 M© pc~^. The average 
value of the virial parameter is a = 5a^R/{GM) ^ 2, 
so that they are gravitationally bound, consistent with 
the observed cloud properties in nucle ar rings of barred 
galaxies (e.g., |McKee &: Ostriker||2007 ). 



5. SUMMARY AND DISCUSSION 

5.1. Summary 

We have presented the results of two-dimensional hy- 
drodynamic simulations of star formation occurring in 
nuclear rings of barred galaxies. We initially consider an 
infinitesimally thin, isothermal gas disk placed under the 
external gravitational potential. The external potential 
consists of a stellar disk, a stellar bulge, a central BH, 
and a non-axisymmetric stellar bar. We do not study 
the effect of spiral arms in the present work. The bar 
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increasing SFR. 



potential is modeled by a Ferrers prolate spheroid with 
the semi-major and minor axes of 5 kpc and 2 kpc, re- 
spectively, and rotates about the galaxy center with a 
patten speed of 33 km s~^ kpc~ . The bar mass, when 
it is fully turned on, is set to 30% of the total stellar 
mass in the spheroidal component, corresponding to a 
strongly barred galaxy. We fix the gas sound speed to 
Cs = 10 km s~^ and the BH mass to 4x 10^ Mq. Our sim- 
ulations incorporate star formation recipes that include a 
density threshold corresponding to the Jeans condition, 
a star formation efficiency, conversion of gas to particles 
representing star clusters or their groups, and delayed 
momentum feedback via SN explosions. To explore var- 
ious situations, we consider both uniform and exponen- 
tial density models, and vary the gas surface density, bar 
growth time, and the total momentum injection in the 
in-plane direction. The main results of this work can be 
summarized as follows. 



The imposed bar potential readily induces a pair of 
dust-lane shocks in the bar region inside the outer- 
most xi-orbit. At early time when the bar potential is 
weak, the dust-lane shocks are placed at the far down- 
stream side from the bar major axis. As the bar poten- 
tial increases, the dust-lane shocks become stronger and 
slowly move toward the bar major axis. The gas passing 
through the shocks loses a significant amount of angu- 
lar momentum, infalls radially along the dust lanes, and 
forms a nuclear ring. The continuous gas inflows provide 
a fuel for star formation in the ring. After the bar po- 
tential reaches its full strength, the dust lanes settle on 
an xi -orbit, while the nuclear ring follows an a;2 -orbit. 
The remaining gas located inside the outermost a;i-orbit 
and outside the dust lanes is gathered to form an elon- 
gated inner ring in the bar region, whose shape is well 
described by an a;i-orbit, as well. Some of the gas located 
outside the outermost a; i -orbit transits to the inner ring 
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Figure 13. Distribution of giant clouds in the nuclear ring of Model U20 at i = 0.36 Gyr. Left: the gas surface density is shown in linear 
scale. Right: cloud boundaries found by the method described in the text are overlaid on the gravitational potential of the gas. 



near the bar ends where xi-orbits crowd. Similarly, the 
gas in the inner ring loses angular momentum when it 
collides with other gas near the bar ends, slowly infalling 
to the nuclear ring through the dust lanes. 

The contact points between the dust lanes and the nu- 
clear ring is a tunnel through which the inflowing gas 
on xi-orbits switches to the a;2-orbit of the nuclear ring. 
About the time when the bar potential is fully turned on, 
the contact points are located near the bar minor axis. 
Over time, the nuclear ring shrinks in size due to the 
addition of low angular momentum gas from outside and 
by collisions of the ring material, which in turn makes 
the contact points rotate slowly in the counterclockwise 
direction. Since the contact points have largest density 
in the nuclear ring, they are preferred sites of star for- 
mation, although star-forming regions can be distributed 
throughout the ring when the mass inflow rate is high. 

The bar potential transports the gas in the bar region 
to the nuclear ring very efficiently, but does not have 
strong influence on the gas orbits outside the bar region. 
This not only makes the bar region evacuated rapidly but 
also reduces the mass inflow rate Mnr dramatically after 
~ 0.3 Gyr. The SFR in nuclear rings displays a single 
primary burst followed by a few secondary bursts before 
becoming reduced to small values. The primary burst 
is associated with the massive gas inflow along the dust 
lanes caused by the growth of the bar potential. The du- 
ration and maximum rate of the primary burst depend on 
the growth time of the bar potential, in such a way that 
a slower bar growth results in a more prolonged and re- 
duced SFR. The secondary bursts are due to the re-entry 
of the ejected and swept-up gas by SN feedback from the 
nuclear ring out to the bar-end regions. Time intervals 
between the secondary bursts are roughly ^ 50 — 80 Myr. 
The peak SFR is attained at i ^ (0.8 — l)Tbarj and the 
duration of active star formation is roughly a half of the 



bar growth time. The SFR is almost equal to Mnr- 
It has a weak dependence on the total gas mass Mnr 
in the ring when SFR ^ 1 M© yr~^, and is not corre- 
lated with AfiMR when SFR ;$, 1 M© yr~^. This suggests 
that star formation in the ring is controlled primarily 
by Mnr rather than Mnr. The relationship between 
the SFR surface density and the surface density of dense 
clumps in nuclear rings found from our numerical models 
are consistent with the usual Kennicutt-Schmidt law for 
circumnuclear starburst galaxies. 

The presence or absence of azimuthal age gradients of 
young star clusters in nuclear rings depends on the SFR 
(or Mnr) in our models. When AInr is small, most of 
the inflowing gas to the nuclear ring is consumed at the 
contact points. In this case, young star clusters that 
form would exhibit a well-defined azimuthal age gradi- 
ent along the ring. When Mnr is large, on the other 
hand, the contact points are unable to transform all of 
the inflowing gas to stars. The extra gas overflows the 
contact points and goes into the nuclear ring. The ring 
becomes massive and forms stars in clumps that become 
dense enough. In this case, no apparent age gradient of 
star clusters is expected since star-forming regions are 
randomly distributed over the whole length of the ring. 
The critical value of A/nr that determines the presence 
or absence of the azimuthal age gradient is estimated to 
be ~ 1 Mq yr~^ in our models, although it depends on 
various parameters such as the ring radius, critical den- 
sity, etc. (eq. [6]). 

Star clusters produced also exhibit a positive radial 
age gradient, such that young clusters are located close 
to the nuclear ring, while old clusters are found away 
from the ring. The primary reason for this is that 
the nuclear ring becomes smaller in size with time, and 
thus star-forming regions gradually move radially in- 
ward. In our models, the radial age gradient amounts 
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to /3 = (ilog(i/yr)/d(r/kpc) ~ 6 — 15. Radial diffusion 
of clusters via mutual gravitational interactions and also 
with the gaseous ring does not affect the radial age gra- 
dient much. 

When the SFR is large (> 1 Mq yr~^), some dense 
clouds are able to grow by accreting surrounding ma- 
terial and form massive clusters. The cluster mass 
function is well described by a power law, with slope 
r = — dlog A''/dlogM* ^2 — 3. A larger slope cor- 
responds to a more massive disk in which more gas is 
available for the ring star formation. When the SFR is 
small, on the other hand, most clusters form near the 
threshold density, leading to a steeper slope with F ^ 3. 
Giant clouds in nuclear rings have typical masses 10^ M0 
and sizes 0.1 kpc. Driven by momentum injection from 
SNe, their one-dimensional internal velocity dispersions 
are supersonic at ^ 20 km s~^. They are gravitationally 
bound with the virial parameter of a ~ 2. 

5.2. Discussion 

We find that the SFR in nuclear rings shows a strong 
primary burst, with its duration and peak value depen- 
dent on the bar growth time, and subsequently a few 
weak and narrow bursts, after which the SFR becomes 
very small. The peak of the primary burst is attained 
roughly when the bar potential is fully turned on. This 
burst behavior of the SFR appears to be a generic fea- 
ture of star formation in nuclear rings of strongly-barred 
galaxies found in numerical simulat ions. For instance, 
A^-body-|- SPH model s prese nted by Heller fc Shlosman 
( [T994| , [fompen et aL| ( pJ95l ), and |Friedn fc^e nz (1995) 
showed that the Hl^'R reaches its peak value, with narrow 
bursts superimposed, about the time when the stellar bar 
fully develops, after which it is reduced to small values. 
In numerical modeling for star forma t ion in the nuclear 
region of the Milky Way, Kim et al. (20111 found that 
the SFR is maximized at i ~ 0.15 Gyr, with a peak value 
~ 0.22 M0 yr~^, and then drops to a relatively constant 
value ~ 0.05 — 0.07 M^ y r~^. The sustained star forma- 
tion inlKim et al. (2011) is thought to arise because the 



Milky Way has a very weak bar that takes a long time 
to clear out gas in the bar region. In this case, the gas 
infall may proceed continuously over an extended period 
of time, a situation similar to the case with a slowly- 
growing bar. 

There is observational evidence that star formation in 
nuclear rings occurs continually over a long period of 
time (a few Gyrs) with successive ^--^ 4 — 1 bursts sep 
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arated by a few tenth s of Gyr s each 

2006| iSarzi et al.|2007[ see also [van der Laan et ai|2013p 

'I'his IS in sharp contrast to our numerical results thai 



sharp contrast to our numerical results that 
show that star formation in nuclear rings is dominated 
by one primary burst before declining to small values, 
with ^0.1 Gyr duration of active star formation. This 
discrepancy is mostly likely due to the fact that our mod- 
els consider only a bar potential for angular momentum 
transport and thus are too simple to describe more com- 
plicated mass inflows in real disk galaxies. If the mass 
inflow rate to a nuclear ring really controls the SFR in 
the ring, as found by our numerical models, the obser- 
vational results call for a need to feed the rings with 
gas episodically for a long time interval. The bar poten- 
tial alone is unlikely responsible for gas supply needed for 
star formation in real nuclear rings. Unless bars are dy- 



namically young, present star formation in nuclear rings 
of nearby barred galaxies requires additional gas feeding. 
One obvious such mechanism is spiral arms that can re- 
move angular momentum at spiral shocks to transport 



gas froni outer disks to the b ar regions (e.g., Lubow et 



al. 



m t 



1986 Kim & Stone 2012), which is not considered 



;he present work. Accretion of halo gas to the disk 



may n ot only rejuvenate bars (e.g.. 



Bournaud & Combes 



2002 1 but also enhance the SFR in the rings (e. g., Jiang &: 
Eimiey|1999{ [Fraternah fc Binney||2006[ |2008[ ). Such gas 
hows niight actually exist, as evidenced by the presence 
of an enhance d number of carbo n stars in the outer spiral 



arms of M33 ( Block et al.ll2007l) . Temporal variations in 
the bar strength (e.g., iBournaud & Combes 200 2J and 
in the bar pattern speed (e..g, [Combes fc Sandersj 1981 1 
are also likely to affect the mass inflow rate to the ring 
and thus the SFR. 



Some galaxies such as IC 4933 (Ryder et al. 2010 ) show 
age gradients of star clusters along the azimuthal direc- 
tion in nuclear rings, while there ar e other galaxies such 
as NGC 7552 dBrandle et^|2012P that d o not show a 



clear age gradient. Mazzuca et al. (2008) analyzed Ha 



data for 22 nuclear rings and found that about half of 
their sample galaxies contain azimuthal age gradients, 
although most of them are not throughout the entire 
ring. They also found that the mean SFR in galaxies 
with azimuthal age gradients is 2.2 ±0.7 M0 yr"-'^, which 
is slightly larger than the mean value of 3. 6 ±1.1 M© yr~^ 
for galaxies with no apparent age gradient. While this 
appears consistent with our numerical results, the large 
dispersions in the mean SFRs suggest that there is no 
fixed SFR that can distinguish between galaxies with 
and without age gradients. In addition, the critical 
SFR for the absence or presence of azimuthal age gra- 
dients is about 1 M0 yr~ ^ in our models , while most 
galaxies in the sample of Mazzuca _gt al. (2008) have 



SFR > 1 Mq yr -^. As equation ^ suggests for the 

maximum SFR, A/^ cPi allowed at the contact points, 
however, the critical SFR depends on many factors that 
may vary from galaxy to galaxy. For example, NGC 
1343 with the most c lear bi-polar age gradient in the 
Mazzuca et al. (2008) sample has the curre nt SFR of 
6.8 M(7) yr~^. Its ring radius is ^ 2 kpc (Comeron 



et al. 2010), which increases the critical SFR by a fac 
tor or 4, assuming that the ring width is proportional to 
the ring size and the other parameters remain the same. 



In addition. 



Af*,CP OC Sr-p oc C? ' 



,3/2 



SO that the level 



of interstellar turbulence and the size of star-forming re- 
gions RsF rnay change M*,cp considerably. 

We find that star clusters that form in nuclear rings 
naturally develop a positive radial gradient of their ages 
owing primarily to the decrease in the ring size over 
^ Gyr in our mod els. This is consistent with the re- 



sults of Jang & Lee (2013) who found that clusters with 
ages ^ 1 Gyr in the nuclear region of NGC 1672 are 
older systematically with increasing radius. Note that 
the radial age gradient holds over a timescale of ^ Gyr 
and may not apply to clusters in a small age r ange since 



the decay o f the ring size is quite slow. Indeed, Mazzuca] 



cay o 



et al. mm found that two (NGC 5953 and 7570) of 



their sample galaxies show a negative radial age gradient 
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of H II region s in the ringsr] Our numerical results plot- 
ted in Figure [iT] also show that when limited to clusters 
with age ~ lO^lO^'^ yr at t = 1 Gyr, younger clusters 
can be found at larger radii, which is due to the stochas- 
tic nature of star formation and ensuing gravitational 
interactions. 

Our results show that the SFR in the nuclear rings is 
tightly correlated with the mass inflow rate to the ring 
rather than the total gas mass in the ring (Fig. [6|. This 
result is seemingly consistent with the results of [Benedict | 
et al. (2002) who found that the SFR in the nuclear ring 
ot a strongly-barred galaxy NGC 4314 is smaller, by a 
factor of 30, than that in a weakly-barred galaxy NGC 
1326 ( Buta et al.||2000 l, even if the gas mass in the ring 
is smaller by only a factor of two. It is interesting to note 
that the gas mass contained in most "gas-rich" nuclear 
rings of barred galaxies in the BIMA SONG sample is 
in a remarkably narrow range of ~ (1 — 6) x 10^ M^) 



(Shethetal, 



2005|P The SFR data presented in Mazzuca 



et al. 1 20U8') combined with the bar strength given m 



Comeron et al. (2010| show that strongly-barred galaxies 
usually have very small present-day SFRs and the SFRs 
in weakly-barred galaxies vary in a wide range, although 
the number of galaxies in their sample is too limited to 
make a conclusive statement. It will be interesting to 
see how the bar strength as well as gas inflows by spiral 
shocks influence the SFR in nuclear rings. 

While star formation is concentrated in nuclear rings in 
our models, observations indicate that star formation in 
some galaxies occurs not only in nuclear rings but also in 
the bar region including dust lanes (e.g., Martin fc Friedli 



19971 IS heth ct al.|[20q0l|Zurita fc Per ez"iiOOH^ 'ijimegreen 
et al.|2 009; Martinez-Garcia & Gonzalcz-Lopezlira 2011). 
While dust lanes themselves are known hostile to star for- 
mation due to strong v elocity shear (e.g. , |Athanassoula| 
1992||Kim et al.|2012a[ ), |Sheth et al.| ( |2000[ ) proposed that 
stars form in inter bar dust spurs in filamentary shape 
that impact t he dust lanes froin the trailing side of th e 



bar (see also Sheth et al. 



20021 [Zurita fc Pe rez"2008p. 



Indeed, plmcgrccn et al. 120U9) interred that some clus 



ters in the nuclear ring of JN GC 1365 actually formed in 
one of the dust lanes by the impact of spurs and subse- 
quently migrated inward to the nuclear ring. The origin 
of these interbar spurs is yet unclear. Apparently, there 
is no filamentary interbar feature in our models. They 
may originate from gas inflows d ue to spir al shocks from 
the region outside the bar (Elmegreen et al.^^2009j , from 
interactions of gas wit h magnetic fields that a re perva- 
sive in the bar region (Beck et al. 1999 2005), and/or 
from other dynamical processes that involve gas cooling, 
self-gravity, etc., which are not considered in the present 
work. It will be an important direction of future work to 
study how spiral arms and magnetic flelds affect the gas 
inflows and star formation in the bar and nuclear regions. 

We gratefully acknowledge I. S. Jang and M. G. Lee for 
sharing their results on the radial age gradient of clus- 
ters found in NGC 1672. We also thank E. C. Ostriker 
and K. Sheth for helpful discussions, and are grateful 

^ Note that NGC 5953 is a non-barred galaxy. 

6 WhiIe |Sheth et al.| ( |2005l l reported that the ring in NG C 6946 
has a mass of ~ 1 0" ^Pi, a higher-resolution observation of |Schin^ 
nerer et al. i 2006 1 gives the ring mass of ^ 4 X 10** Mq . 
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